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Spin Caloritronics in graphene with Mn 
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We show that a device made of graphene with Mn adatoms trapped at single vacancies feature spin-dependent 
Seebeck effect, thus enabling the use of this material for spin caloritronics applications. More than that, it 
is possible to tunc its thermoelectric properties with a gate potential in a way the system presents either 
a nearly total spin polarized current, flowing in one given direction, or currents for both spins flowing in 
opposite directions without net charge transport. 

PACS numbers: 72.20.Pa, 72.80.Vp, 85.75.-d, 72.25.Ba 
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The field of spin caloritronics deals with the interaction 
of spin and heat currents, that is, the coupling between 
spintronics and thermoelectrics 1 . Graphene 2,3 is a po- 
tential candidate material for spintronic devices 4 ' 5 due 
to its long mean free path and weak spin-orbit coupling. 
Furthermore, it has been shown that graphene has poten- 
tial for thermoelectric applications, both theoretically 6 ' 7 
and experimentally 8-10 . Therefore, at least in princi- 
ple, graphene is a good candidate for real applications 
in spin caloritronics. Devices made of graphene nanorib- 
bons with zigzag edges have already been proposed to 
this end 11 ' 12 . However, it is an experimentally difficult 
task to synthesize nanoribbons with well defined zigzag 
edges 13 . Bi-dimensional graphene itself, although a good 
conductor for spin polarized currents, is spin degener- 
ate, so it cannot be used as a source for spin polarized 
currents. One way to lift its spin degeneracy that has 
been investigated recently is to dope graphene with metal 
adatoms 14-19 . Although their tendency is to diffuse and 
to form clusters, they can be trapped in defects like single 
vacancies (SV), where they are highly stable 15 ' 20 . 

It has already been shown by some of us that Mn atoms 
trapped at SV (Mn@SV) in graphene are tightly bound 
(with a binding energy of —6.0 eV) and act as a source 
of spin polarized currents under a bias voltage. More- 
over, these currents can be tuned by a gate potential 15 , 
V g . Being that so, it is natural to ask if the thermoelec- 
tric properties of this kind of device would also be spin 
dependent. In this work, wc address this question by cal- 
culating the thermoelectric properties and currents of a 
device consisting of graphene with a Mn@SV under the 
effect of a gate voltage and with contacts held at different 
temperatures (see Fig. 1). 

Our calculations have been performed with non- 
equilibrium Green functions coupled to density func- 
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FIG. 1. (Color online) Schematic view of the geometry of the 
device studied, a graphene sheet (periodic in the x direction) 
with a Mn atom trapped at a single vacancy, coupled to pris- 
tine graphene contacts. The central region of the device is 
under the effect of gate potential. 



tional theory (NEGF+DFT) in the linear regime. First, 
the geometries, comprised of a Mn@SV in a graphene su- 
percell (13.0 x 20.0 x 32.5 A 3 ), have been fully relaxed 
with a force criterion of 0.02 eV/A using the siesta 
code 21 . Then, we calculated the transmittances using the 
TRANS AMPA2 code 22 . The electrodes were considered to 
be semi-infinite pristine graphene sheets. In all calcula- 
tions we used the Perdcw-Burke-Ernzerhof generalized 
gradient approximation 23 for the exchange-correlation 
functional, norm-conserving pseudopotentials 24 and a 
double- £ polarized basis. We used an energy mesh cut- 
off of 300 Ry and a 7 x 1 x 4 fc-point sampling, in 
the Monkhorst-Pack 25 scheme, to integrate the Brillouin 
zone. For the electronic transport calculations we em- 
ployed 101 fcj^-points (x direction in Fig. 1). The gate 
potential was simulated by adding a smooth electrostatic 
potential to the Hamiltonian in a finite region (a xy slab) 
containing the Mn@SV in the self-consistent cycle, thus 
allowing screening effects by charge rearrangement 15 . 
This potential consisted of a plateau, U(z) = eV g for 
zq < z < Z\, represented by a dashed line in Fig. 1, with 
smooth sinusoidal steps with width t, indicated by the 
dotted lines in Fig. 1. U(z) = elsewhere. In the calcu- 
lations we defined the gated region trough zq = 11.24 A, 



2 




V (V) 

g 

FIG. 2. (Color online) Seebeck coefficients for the up and 
down spin channels versus gate potential for an average tem- 
perature of 300 K. The lines are just guides to the eyes. 

zi = 21.23 A and t = 1.25 A. A vacuum layer of 20 A 
was used to avoid spurious interactions between periodi- 
cal images in the direction perpendicular to the graphene 
plane (y in Fig. 1). 

In the Landauer-Buttiker formalism, the current is 
given by 26 

h = \ J %(E; V g ) (f L - f R ) (E - //) dE, (1) 

where a I is the spin, e is the electron charge, h is the 
Planck constant, T„{E\ V g ) is the spin resolved transmit- 
tance function, which depends on the gate voltage V g \ 
fum = f(E,n,T L ( R ^) is the Fermi-Dirac distribution 
function, /i is the chemical potential of the electrodes 
and Tl(r) is the temperature of the left(right) lead. 

The Seebeck coefficient, also named thermoelectric 
power, is a measure of the voltage induced by a tem- 
perature difference on a device or material and is defined 
as S = - A V/ AT 1 1 =0 . In the limiting case of AV 
and AT — > (the linear regime) an expression for S can 
be derived from an expansion of Eq. (1), given by 27 

1 J%(E;V g )(-§i)(E-^dE 
& fTAE;V g )(-§l)dE 

where T is the average temperature between the contacts. 

In Fig. 2 we show the calculated Seebeck coefficients 
for the up and down spin channels as a function of V g 
at room temperature (T = T L (r\ = 300 K). Close to 
zero gate, around V g = 0.24 V, where Stot (= St + Sj.) 
crosses V/K; the system presents St and Sf with ap- 
proximately the same absolute value (~ 12 //V/K) and 
opposite signs, meaning that the thermally induced volt- 
ages generated at a given temperature difference will 
drive the electrons with opposite spins in contrary di- 
rections. In the range V g = 0.1 — 0.4 V, St shows a 
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FIG. 3. (Color online) Current versus gate voltage, V g , for a 
fixed temperature difference of AT = 30 K (the left electrode 
is kept at room temperature, Tl = 300 K, and the right 
electrode at Tr = 330 K). The lines are just guides to the 
eyes. 



nearly constant behavior, decreasing very slowly, from 
-14.2 nV/K at V g = 1 V to -15.6 /xV/K at V g = 4 V. 
Conversely, S^ shows a parabolic shape, varying around 
zero with small magnitude values. Thus, in this range 
of gate voltages, the system is in a regime presenting 
|S t | ~ 14 - 15 fiV/K and Sj. - 0. On the other hand, as 
V g decreases to negative values, St grows steadily from a 
negative value of —11.8 //V/K at V g = V to a positive 
one, +3.6 //V/K at V g = — 4 V, crossing zero at a gate 
voltage close to V 9 = -3 V; Sj, grows faster and shows 
a peak around V g = —2.5 V, where Sj, = 83.2 //V/K and 
St = —1.9 //V/K. Then, close to the peak maximum, 
where Stot and Sj, intersect, the system presents high S^ 
and negligible St- Therefore, we identify three distinct 
regimes into which the present device can be tuned: (i) 
negligible St and high positive Sj, (at negative V g ); (ii) 
St ~ — Sj_, close to zero V g ; and (iii) negative St and 
negligible S^ (at positive V g ). 

To further characterize our device, we calculated the 
current (Eq. 1) as a function of the gate voltage for 
a fixed temperature difference of 30 K between the left 
and right electrodes (fixing Tl = 300 K and Tr = 330 
K) while keeping the external bias voltage at zero. The 
results, plotted in Fig. 3, show more clearly how the 
current varies as V g changes, passing through the three 
distinct regimes discussed. For the current, which can 
be approximated by I a = ~G a S a AT ', where G a is the 
conductance, for small AT; the values or ranges of gate 
voltages at which every regime happens can be somewhat 
different when comparing to S due to the dependence of 
the current on the conductance as well. However, the 
behavior of the currents resembles that of the Seebeck 
coefficient. Regime (i) (If ~ and high Tj,) occurs in 
a narrower range around V g — —2.5V due to the faster 
increase of I\_ with decreasing V g . Regime (iii) (negative 
If and If ~ 0) happens in the range V g = 2 — 4 V and 
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FIG. 4. (Color online) Current versus temperature difference 
(Tl = 300 K and T R = T L + AT) for the applied gate voltages 
of V g = -2.5 V, 0.4 V and 2.0 V, from (a) to (c) respectively. 
Notice the different scales for the current. The lines are just 
guides to the eyes. 



regime (ii), at which there is no net charge transport 
(If = —If), at V g = 0.4 V, presenting a slight shift from 
the value of V g = 0.24 V for the Seebcck coefficient. 

The behavior of the current at different temperature 
differences for the gate voltages V g = —2.5, 0.4 and 2.0V, 
which are representative of each regime described above, 
(i) to (hi) respectively, is shown in Fig. 4 for a fixed 
Tl = 300 K and varying the temperature of the right 
contact as Tr = Tl + AT. As can be seen in the figure, 
the currents follow a linear trend with AT, indicating 
that the system is in the linear transport regime. For 
V g = —2.5 V, when the device is in regime (i) (see Fig. 
4a) , If w A/m in all the calculated range of tempera- 
tures whereas I± reaches up to 1.69 A/m at AT = 40 K. 



Then in this case, the total current, L. 
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sists practically of down-spin electrons only, and has a 
charge polarization, given by (\If\ - \If\)/{\If\ + \If\), 
of —98%. At V g = 0.4 V the system is in regime (ii) 
(Fig. 4b), where for all the calculated ATs, the currents 
for both spins are similar in magnitude (~ 0.52 A/m 
at AT = 40 K) but have opposite directions (given by 
the sign of the currents). The up- and down-spin cur- 
rents cancel each other out and Itot is then negligible. 



Thus, the device displays counter-propagating spin cur- 
rents without any net charge transport. At V g = 2.0 V 
the system reaches regime (iii), showing a behavior op- 
posite to (i). The down-spin current is almost negligible, 
Jj. ks 0, in all the temperature range calculated and If 
attains -0.50 A/m for AT = 40 K. The total current is 
then basically made of up-electrons only, 7 tot « If, and 
shows a charge polarization as high as 94%. 

It is important to note that our proposed device is quite 
small and that a more realistic one would have a larger 
area with Mn randomly distributed over it, and possi- 
bly trapped at different types of defects other than only 
single vacancies. However, if the defects are created by 
ion bombardment, for example, the predominance of the 
type of defects can be adjusted by the beam energy 28,29 . 
Also, anti-ferromagnctic coupling among neighboring Mn 
should be prevented. This could be avoided either by 
keeping a low enough density of Mn@SV or by aligning 
the Mn local magnetic moments using a magnetic field 
(or both). Under these conditions, the effect of disor- 
der in the distribution of defects may lead the system 
into a localized regime at which the spin asymmetry of 
transport properties, the Sccbeck effect included, could 
be enhanced 30,31 . Other point worth considering is the 
use of other metals as dopants (such as Fe, Ni, Co or Au 
and Ag). In principle, any metal showing spin-dependent 
transmittance close to the Fermi level may present spin- 
dependent Seebeck effect (see Eq. 2). Finally, the effects 
of the presence of different substrates should be taken 
into account when building the real device. Depending 
on the substrate, graphene may be p- or n-type doped 
already and then, the gate voltages needed to attain the 
different regimes will be somewhat different. Yet, the 
overall picture for the device will still be the same. Local 
substrate doping can also be used as a substitute for the 
effect of the gate potential, raising or lowering the energy 
levels of the region containing the Mn atoms. 

In summary, we have proposed a proof-of-concept de- 
vice made of graphene with Mn atoms trapped in single 
vacancies showing a highly tunable spin dependent See- 
beck effect, which could open a new door for applications 
of graphene in spin caloritronics. 



ACKNOWLEDGMENTS 

We would like to thank A. R. Rocha for a critical read- 
ing of the manuscript and useful suggestions. This work 
has received financial support from the brazilian agencies 
CNPq, INCT-Nanomateriais de Carbono and FAPESP. 
The calculations have been performed in the GPU cluster 
at CENAPAD-SP. 

1 G. E. W. Bauer, E. Saitoh, and B. J. van Wees, Nature Mater. 
11, 391 (2012). 

2 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, 
S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 306, 
666 (2004). 

3 A. H. Castro Ncto, F. Guinea, N. M. R. Peres, K. S. Novoselov, 
and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 



4 



4 V. M. Karpan, G. Giovannetti, P. A. Khomyakov, M. Talanana, 

A. A. Starikov, M. Zwierzycki, J. van den Brink, G. Brocks, and 

P. J. Kelly, Phys. Rev. Lett. 99, 176602 (2007). 
5 N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and B. J. 

van Wees, Nature 448, 571 (2007). 
6 D. Dragoman and M. Dragoman, Appl. Phys. Lett. 91, 203116 

(2007). 

7 B. Zhou, B. Zhou, Z. Liu, and G. Zhou, J. Appl. Phys. 112, 
073712 (2012). 

8 Y. M. Zuev, W. Chang, and P. Kim, Phys. Rev. Lett. 102, 

096807 (2009). 
9 D. Basko, Science 334, 610 (2011). 

10 X. Wu, Y. Hu, M. Ruan, N. K. Madiomanana, C. Berger, and 

W. A. de Hecr, Appl. Phys. Lett. 99, 133102 (2011). 
n M. Zeng, Y. Feng, and G. Liang, Nano Lett. 11, 1369 (2011). 
12 M. Zeng, Y. Feng, and G. Liang, Applied Physics Letters 99, 

123114 (2011). 
13 F. Schwierz, Nature Nanotech. 5, 487 (2010). 
14 0. Cretu, A. V. Krashcninnikov, J. A. Rodriguez-Manzo, L. Sun, 

R. M. Nieminen, and F. Banhart, Phys. Rev. Lett. 105, 196102 

(2010). 

15 M. P. Lima, A. J. R. da Silva, and A. Fazzio, Phys. Rev. B 84, 
245411 (2011). 

16 E. J. G. Santos, A. Ayuela, and D. Sanchez-Portal, New J. Phys. 

12, 053012 (2010). 
17 C. Cao, M. Wu, J. Jiang, and H.-P. Cheng, Phys. Rev. B 81, 

205424 (2010). 

18 E. J. G. Santos, D. Sanchez-Portal, and A. Ayuela, Phys. Rev. 



B 81, 125433 (2010). 
19 A. V. Krashcninnikov, P. O. Lchtinen, A. S. Foster, P. Pyykko, 

and R. M. Nieminen, Phys. Rev. Lett. 102, 126807 (2009). 
20 Y. Tang, Z. Yang, and X. Dai, J. Chem. Phys. 135, 224704 

(2011). 

21 E. Artacho, D. Sanchez- Portal, P. Ordejon, A. Garcia, and J. M. 
Soler, Phys. Status Solidi (b) 215, 809 (1999). 

22 F. D. Novacs, A. J. R. da Silva, and A. Fazzio, Braz. J. Phys. 36, 
799 (2006); J. E. Padilha, M. P. Lima, A. J. R. da Silva, and 
A. Fazzio, Phys. Rev. B 84, 113412 (2011); The code has been 
modified to allow the use of graphical processing units (GPUs). 

23 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 
3865 (1996). 

24 N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 (1991). 

25 H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976). 

26 R. Landauer, Philos. Mag. 21, 863 (1970); M. Biittiker, Y. Imry, 
R. Landauer, and S. Pinhas, Phys. Rev. B 31, 6207 (1985); 
M. Biittiker, Phys. Rev. Lett. 57, 1761 (1986). 

27 U. Sivan and Y. Imry, Phys. Rev. B 33, 551 (1986). 

28 A. W. Robertson, B. Montanari, K. He, J. Kim, C. S. Allen, 
Y. A. Wu, J. Olivier, J. Neethling, N. Harrison, A. I. Kirkland, 
and J. H. Warner, Nano Letters 13, 1468 (2013). 

29 A. W. Robertson, C. S. Allen, Y. A. Wu, K. He, J. Olivier, 
J. Neethling, A. I. Kirkland, and J. H. Warner, Nature Commu- 
nications 3 (2012), 10.1038/ncomms2141. 

30 A. R. Rocha, T. B. Martins, A. Fazzio, and A. J. R. da Silva, 
Nanotechnology 21, 345202 (2010). 

31 J. M. de Almeida, A. R. Rocha, A. J. R. da Silva, and A. Fazzio, 
Phys. Rev. B 84, 085412 (2011). 



